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Abstract 

We present a theory of hysteretic phenomena in Bose gases, using superfluidity in one dimensional 
rings and in optical lattices as primary examples. Through this study we are able to give a 
physical interpretation of swallowtail loops recently found by many authors in the mean-field 
energy structure of trapped atomic gases. These loops are a generic sign of hysteresis, and in the 
present context are an indication of superfluidity. We have also calculated the rate of decay of 
metastable current carrying states due to quantum fluctuations. 
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I. INTRODUCTION 



Quantum degenerate bosonic atoms have proven important for studying macroscopic 
quantum phenomena (for a review see [0]). The order parameter of the condensed phase 
is a macroscopic quantum wavefunction which, unhke single particle wavefunctions, can be 
directly probed in an experiment. The interplay between this macroscopic wavefunction 
and interactions leads to a variety of novel effects, the most well known of which is su- 
perfluidity. Here we explore superfiuid phenomena in a dilute atomic gas with short range 
interactions. As we will show, superfluidity is naturally viewed as a hysteretic response to 
rotation, motivating a more general study of hysteresis. 

In our study of superfluidity, we quantify the roles played by interactions, finite size 
effects, and impurities in the behavior of a weakly interacting gas of one dimensional Bosons, 
showing that persistent currents can exist when the interactions are strong compared to any 
impurity potentials, but weak enough to not produce large phase fluctuations. We present a 
detailed discussion of the energy landscape of such gases, reveahng a non-trivial topography. 
In the limit of weak interactions we calculate the lifetimes of persistent currents. 

In addition to gaining insights into superfluidity within a one dimensional geometry, our 
approach provides an intuitive understanding of swallowtail energy loops found in mean-field 
studies of Bose gases within periodic potentials 0, |], |[. We show that such loops are a 
generic feature of hysteresis and, in the case of atoms in a periodic potential, the loops 
are a manifestation of superfluidity. We discuss the underlying quantum scaffolding which 
supports this mean-field structure, and identify other settings where it can be observed. 

In section I A we introduce the basic phenomenon of hysteresis. The remainder of the 
introduction discusses superfluidity and provides examples of scenarios in which a Bose con- 
densate will behave hysteretically. Section |I| discusses microscopic models for superfluidity 
in both a ring shaped geometry and in an optical lattice. The remainder of this paper 
analyzes these models. 

A. Generic properties of hysteresis 

Hysteresis is the phenomenon where the state of a physical system depends upon its 
history. The canonical example is a ferromagnet, which in zero applied magnetic field has a 
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spontaneous magnetization, conventionally taken to be in the z direction. This magnetiza- 
tion is robust in that it is not significantly changed applying a small field in the —z direction. 
However, if a strong enough field is applied, the magnetization can 'fiip', and point in the —z 
direction. When the applied field is reduced to zero, the magnetization does not revert to 
its original orientation, but remains pointing in the —z direction. In this example, and the 
ones that follow, we see that the response of the system lags behind the applied stimulus. 

In a classical system, hysteresis is conveniently thought about in analogy to the Landau 
theory of phase transitions [|5|. One considers the property of interest (in this case the 
magnetization M) to be an order parameter. An energy landscape is produced by calculating 
the energy of the system as a function of this order parameter. The applied field (here the 
magnetic field if) changes this landscape. 

Hysteresis occurs when the energy landscape has more than one minimum, as depicted in 
FIG. |l|a (for similar figures calculated within a microscopic model, see FIGS. |12| and |13]) . For 
example, both magnetization in the z and —z directions might be local minima. Applying a 
field tilts the landscape, and reduces the barrier. At some critical field, the barrier disappears 
and the system jumps into the global minimum (FIG. |I]b and c). The phenomenon where 
the barrier disappears goes under several names; in the theories of phase transitions and 
of gradient dynamics |]^ it is respectively known as a spinodal or a catastrophe. In more 
mathematical treatments it is referred to as a "saddle-node bifurcation" . 

Figure ^ gives a generic depiction of the energy of the extrema of the energy landscape 
(again, similar figures calculated from microscopic models are shown in FIG. |ll]). A dis- 
tinctive loop is seen. This loop, referred to as a "swallowtail" by Diakonov et al. is a 
general feature of the spectrum of a hysteretic system. It exists because for some range of 
fields there are three extrema (two local minima and a maximum). At the point labeled by 
(c), one of the local minima meets up with the maximum, and they annihilate one another. 

To better match the dynamical systems literature, it would be preferable to not refer to 
FIG. |]as a swallowtail, and instead reserve the term for the similar structure in FIG. Q that 
occurs when one has two control parameters. The second control parameter changes the size 
of the loop, and can be tuned so that the loop, and all hysteresis, vanishes. A model which 
gives rise to this latter structure will be discussed later. 

A possible point of confusion here is that the term "swallowtail" is traditionally used to 
discuss not the energy structure, but rather the catastrophe set, which is the points (in the 
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control parameter space) where the number of extrema in the energy landscape change. The 
catastrophe set corresponding to FIGS. |^ and ^ respectively consists of two points and the 
cusp-like structure in FIG. |[ Thus the swallow-tail energy spectrum is associated with a 
cusp catastrophe and not a swallow-tail catastrophe. 

The local minima in the energy landscape are of great physical importance, as the system 
typically resides in their vicinity. Saddle points, and local maxima, are also important in 
that the rate of transitions from one minima to another are governed by the lowest barrier 
separating the minima. In classical systems, these transitions are typically caused by thermal 
fluctuations, and occur at a rate proportional to e'^^^'^'''^, where Eb is the barrier height, kb 
is Boltzmann's constant, and T is the temperature. It should be noted that only in very 
rare physical situations does the system spend much time at one of these extrema. 

In a quantum mechanical system the scenario for hysteresis that we have discussed be- 
comes more complicated. The basic difficulty is that the order parameter is generally not a 
constant of motion. In this case one does not know how to answer questions like "what is 
the energy of the system when the magnetization points in the +i direction?" There may 
simply not exist any energy eigenstates for which the magnetization points in that direction. 
Consequently it is by no means obvious how to construct an energy landscape, and what 
significance it will have. 

There are three, roughly equivalent, methods of circumventing this difficulty. The first 
approach is to write the Hamiltonian as a sum of two terms, H = H^iag + H', where the order 
parameter commutes with i^diag- This diagonal term is the projection of the Hamiltonian 
into the space where the order parameter has a definite value. For example, if we have a 
spin system where the z component of the magnetization is our order parameter, then -ffdiag 
would be diagonal in a basis US', 6*2)}, where 5* is the total spin, and Sz is the projection 
of the spin along the z axis. If H' is small, one can neglect it for the sake of drawing the 
energy landscape. The second approach is to use a variational scheme, where one writes 
down "reasonable" wavefunctions which are parameterized by the order parameter. The 
expectation value of the energy in these states is an approximation to the energy landscape. 
The final approach is to use a mean-field theory in which the order parameter is a constant 
of motion. This discussion will be more concrete once microscopic models are introduced in 
section ^ and used to produce energy landscapes. 

All three of these schemes share the feature that if the system starts in a local minimum of 
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the energy landscape, there can be matrix elements in the original Hamiltonian which allow 
the system to tunnel to another minimum. This procedure can be thought of in analogy to 
classical thermally activated transport, where due to thermal fluctuations the system can 
jump from one minimum to another. Here it is quantum fluctuations which allow the system 
to move between minima. 

B. Superfluidity 

We now turn to a discussion of superfluidity, a phenomenon which manifests itself in many 
related ways, including, dissipationless flow, quantized vortices, reductions in the moment 
of inertia, and the existence of persistent currents. We focus on the latter phenomenon, 
which was flrst observed in ^He [0. In an idealized version of these experiments, an annular 
container of helium is rotated while cooling to below the lambda temperature, where it 
becomes a superfluid. When the container is then stopped, one observes that the fluid 
continues to rotate - maintaining its velocity for extremely long times. Arguments based 
solely on Galilean invariance show that this current carrying state cannot be the ground 
state of the system, and is therefore an extremely long lived metastable excited state ^j. It 
is observed that the lifetime of these currents decrease with increasing velocity, and there is 
a critical velocity fc, above which no persistent currents exist. 

For our purposes it is convenient to think of such currents in terms of a hysteretic response 
to rotation. Imagine starting with an annular container of superfluid which is at rest. If the 
container is slowly rotated in a clockwise direction the fluid remains at rest (in the rotating 
frame this is a persistent current). If one rotates faster and faster, the relative velocity 
between the container and the fluid eventually exceeds the critical velocity, excitations are 
formed, and the fluid accelerates. At this point a persistent current has developed in that 
even if one stops rotating the container then the fluid will continue to flow. The flow can be 
stopped if one rotates the container sufficiently fast in a counter-clockwise direction. The 
principle is simply that when the relative velocity between the fluid and the container exceeds 
Vci the fluid accelerates. Thus the fluid flow lags behind the applied rotation, resulting in 
the hysteresis loop sketched in FIG. ^. By the arguments of section I A, one must therefore 
see energy structures analogous to those in FIG. |l] and||. 

Here we wish to understand the origin of this dramatic effect from a microscopic model. 
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Standard descriptions of superfluidity ^ attribute the long life of these currents to the 
scarcity of low energy excitations. In the present setting, it is more natural to think of 
superfluidity in terms of the ability of the fluid to screen out impurities. The basic argument, 
which will be given in more detail later, is that an ordinary (non-super) fluid does not 
support persistent currents because the fluid particles scatter off of small imperfections in 
the walls of the container, exchanging angular momentum and eventually equilibrating with 
those walls. Through a collective deformation of the macroscopic wavefunction, quantum 
degenerate bosons are able to screen out the imperfections. Since the fluid effectively see 
smooth walls, it does not slow down. In the body of this paper these imperfections will be 
modeled as an 'impurity' potential. 

From this physical picture, one can anticipate many of our results. In particular, there 
are two natural control parameters, the rate of rotation and the strength of interactions 
relative to the impurity potential. We will flnd an energy structure similar to FIG. where 
these two control parameters correspond to H2 and Hi. 



C. Optical lattices 

Superfluidity is not limited to a ring geometry. As we explain, superfluid properties nat- 
urally appear for Bose particles within a periodic potential. Due to their importance in 
solid-state physics, quantum phenomena in periodic potentials are very well studied theo- 
retically and there has been rapid progress on experimental studies of Bosons in periodic 
potentials, where the periodicity is produced using standing waves of light (optical lattices) 
(for a review see |]10|)- Many of the single particle phenomena of solid-state physics have 



been observed in these artiflcial lattices, including band structure, Bloch oscillations, and 
Zener tunneling |T^. These solid state concepts are reviewed below, and play an important 
role in our discussion of superfluidity. Further theoretical discussions of these phenomena 
in cold gases and relevant discussion of how interactions screen the lattice can be found in 



TT| . Although not directly related to our study of hysteresis, it is worth mentioning that 
correlated many-body states, such as Mott insulators |]l2l, have been observed in atoms 
trapped in an optical lattice. 

Here we use superfluidity to reexamine theoretical studies of mean-fleld energy loops 
in of atoms in optical lattices 0. We understand the key features of these studies by 
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starting with the energy structure of the non-interacting single particle states. As discussed 
in textbooks |13[, the states available to noninteracting particles in a periodic potential 
are labeled by two quantum numbers, a band index u, and a crystal momentum k. The 
wavefunctions of these states are of the Bloch form 

Mr) = e^'-'-v^ir), (1) 

were v,y{r) shares the periodicity of the lattice, and k is restricted to the first Brillouin 
zone. Limiting our discussion to one dimension with lattice periodicity L, the first Brillouin 
zone corresponds to momenta \k\ < tt/L. For simplicity we use dimensionless units where 
L = 1. In FIG. ^, the lowest energy band is sketched in an extended zone scheme, where the 
energy is extended periodically to k's outside of the first Brillouin zone. This periodicity is 
the source of the phenomenon known as Bloch oscillations. Imagine starting with a single 
particle in the k = state. If an external force is applied to the particle, it will accelerate 
and k will increase. For sufficiently weak acceleration, the state will adiabatically follow the 
solid curve in FIG. |^. When k has increased to 27r, the system has returned to its initial 
state. Thus a constant force leads to periodic oscillations. If the force is too strong, the 
adiabaticity condition is violated, transitions are made to higher bands, and one no longer 
sees the Bloch oscillations. This breakdown is known as Zener tunneling. 

A similar scenario can be considered for Bose condensed atoms. In the ground state 
all of the particles reside in the lowest energy Bloch state. Like the single particle case, 
when a force is applied, the crystal momentum k increases. However, as a superfluid, 
the condensate is able to screen out the periodic potential. Thus, for sufficiently strong 
interactions, instead of following the solid curve in FIG. the system follows a path closer 
to the dashed curve, corresponding to the spectrum of states in the absence of the periodic 



potential. The microscopic model which will be introduced in |IIB| confirms this picture, and 
one can identify the set of states visited during this adiabatic acceleration as local minima 
in a mean-field energy landscape. When the fiuid's velocity exceeds the critical velocity, it 
looses the ability to screen the lattice. Thus the energy curves terminate at some point, and 
energy extrema take on the structure in FIG. |^, where one has a crossing of local minima. 
One minimum corresponds to the fiuid moving to the right, the other to fiuid moving to the 
left. These two states have different momentum, but share the same crystal momentum. For 
purely topological reasons, the presence of two local minima at a given value of k guarantees 
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that there is a saddle point separating them. This barrier state is also shown in FIG. ^ as 
the dotted line forming the 'top of the swallow's tail'. As will be discussed in SEC. |111 B| , 
the barrier state corresponds to a 'phase slip'. 

The question of what happens to the system when the cloud is accelerated past the 'end 
of the loop' (marked by a grey circle on the figure) is discussed by Wu et al. 0, |T^. Clearly 
adiabaticity must break down at this point, and crossing this point from left to right, then 
back again will not return the system to its original state. Thus some sort of hysteresis has 
developed. 

Note that as the interaction strength is reduced the loops in FIG. |^ become smaller and 
eventually disappear. Thus if one identifies Hi with k and H2 with the interaction strength, 
the mean-field energy extrema near k = n must have the full swallowtail structure shown in 

FIG. a 



D. Josephson junctions 

We conclude the introduction by discussing a hysteretic Bose system in which the hys- 
teresis is not associated with persistent currents, namely, a gas of particles with attractive 
interactions in a double- well trap as depicted in FIG. The control parameters here are 
the strength of interactions and the bias 6 which is applied between the two wells. 

FIG. P illustrates the transformations which give rise to a hysteresis loop in this system. 
One starts with the left well of much lower energy than the right {6 > 0). The ground state 
consists of all of the particles bunched up on the left. The bias is then slowly decreased, 
and made slightly negative, so that the right hand well has lower energy. In the true ground 
state all of the particles are sitting in the right hand well. Nonetheless, the particles actually 
stay in the left hand well. This behavior is understood by noting that in order to move the 
particles from the left hand well to the right one has to first move a single particle. Although 
such a move saves the potential energy of the bias, separating that one particle from the 
others makes the interaction energy less negative. For small enough bias moving a single 
particle increases the total energy, and the state with all of the particles sitting in the left 
hand well is a local minimum of the energy. If the bias is made more negative, the potential 
energy savings of moving a particle to the right hand side eventually becomes greater than 
the interaction energy cost. The particles then all jump to the right hand well. The whole 
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process can be reversed, and a hysteresis loop is formed. 

For weaker interactions, the value of 6 at which the metastable state disappears becomes 
smaller. For sufficiently weak interactions, no hysteresis occurs. Thus identifying H2 with 
the interaction strength, and Hi with the detuning 6, the energy landscape has extrema 
matching FIG. |^. 

Experimentally such a double-well trap can be formed by considering two cells in an 
optical lattice, or by carefully arranging magnetic and optical fields as in PI. We will not 
explicitly discuss models for this system, as there exist many excellent treatments in the 
literature [llH 



II. MICROSCOPIC MODELS 



We now construct microscopic models of the superfiuid systems described in the introduc- 
tion. In addition to verifying the qualitative structures already discussed, these models allow 
us to make quantitative predictions about the behavior of a gas of bosons. In particular, as 
mentioned in section I A, a quantum mechanical system can tunnel from one minimum in 
the energy landscape to another (which would, for example, lead to the decay of persistent 
currents). We are able to calculate the rate of such tunneling. 



A. Superfluidity 

As a microscopic model of persistent currents, we study a one dimensional ring of length 
L, rotating at frequency Q, containing a cloud of bosons of mass m which interact via short 
range interactions. Measuring energy in terms of /2mL'^, {% is Plank's constant) the 
Hamiltonian in the rotating frame is 

H = ^(27rj + ^fc]c, + {g/2) ^ cJ4qc^, (2) 

j j+k=l+m 

where $ = 2mL'^Q/h and g > are dimensionless measures of the rotation speed and the 
strength of the interactions. Operators c], create bosons with angular momentum jh. The 
model @ could be experimentally realized by cooling an atomic gas in an annular trap with 
harmonic confinement of frequency uj± to such an extent that only the lowest transverse 



mode is occupied (for recent experimental progress on annular traps see IT^)- As long 
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as the trap length d± = yfi/muj± is larger than the scattering length the interaction 
parameter would then be given hj g = ATiagL/ d\. 

Despite its apparent simplicity, this one dimensional model is quite rich. It is a canonical 
example of a Luttinger liquid [|I^ whose behavior can be studied via the Bethe ansatz [|19| . 
Two properties worth noting are: 1) at (7 = it describes a non-interacting Bose gas; and 2) 
at (7 00 it can be mapped onto a gas of non-interacting fermions. In neither of these limits 
is the system superfluid, however we show that for small positive g the system is superfluid. 
Here we will study how this superfluidity develops as g is tuned from 0, finding the structure 
discussed in section IB. The equally interesting question of how this superfluidity breaks 
down as (7 — > cxD will not be discussed. 

The model Hamiltonian (Q) is invariant under rotation, and therefore conserves angular 
momentum. A trivial consequence is that if a current is started in this system it will never 



decay. Thus, as aptly pointed out by Kagan et al. to study superfluidity one must 

add an impurity which breaks the symmetry. In an experimental setting such terms are 
always present due to imperfections in the apparatus. It is quite instructive to imagine 
artificially introducing such an impurity (for instance by using a laser which interacts with 
the atoms via dipole forces), and being able to control its strength. Conventional discussions 
of superfluidity focus on ^He, which is strongly interacting, and whose behavior is largely 
insensitive to the strength of the impurities. In a weakly interacting setting (especially in 
1-D) this is no longer the case, and the strength of the symmetry breaking term is extremely 
important. The system's behavior is relatively insensitive to the exact form of the impurity. 
Two natural models are a point scatterer H^^t = XJ2kqclcq and a sinusoidal potential 
Hsin = XJ2k (cfcCfc-i + clck+ij- In both cases A measures the strength of the perturbation. 



B. Optical lattices 

A model for particles in a periodic potential can be constructed which has the same 
structure as (0) with an impurity. The rotation speed $ and the impurity potential are 
respectively mapped onto the crystal momentum and the lattice potential. 

Introducing the field operator '?/'(x), which annihilates a particle at position x, the Hamil- 
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tonian for particles in one dimension interacting with a local interaction is 

H = fdx ^'^^^ • + V(xU^ij + ^tZ-VV^, (3) 
J 2m 2 

where V{x) is the lattice potential, g parameterizes the interactions, and the argument 

X is assumed for each of the field operators. The periodic potential can be written as 

V{x) = J2j e'^P^^Vj, where pj = jki = 27rj/L are reciprocal lattice vectors. 

In analogy to (|I]), it is convenient to write our field operators in a Bloch form, 

"^'^^^^ (4) 

where ipj{k) is the Bose operator which annihilates a particle with momentum h{k + 2Txj) / L. 
These obey the standard commutation relationships \il)j{k),ip^y{k')] = 6jj'6{k — k'). In this 
decomposition, k is the (dimensionless) crystal momentum which runs from — tt to vr and pj 
are reciprocal lattice vectors. The index j plays the role of a band index in the limit that 
the lattice potential vanishes. 

In terms of these operators, the Hamiltonian takes the form 

^ = ^ (^H„ + i/pot + HlT^ + H[^:) (5) 

i^kin = /rfA;E,(2vrj + fc)VM (6) 

i^pot = /t^A;E,,KV'J+,V'- (7) 

Hl7'^ = Jdk {g/2) Y.n+n=n+H AA^n^u (») 

HLt = Id{k}ig/2)E{n44^3^^ (9) 

where the respective terms in represent kinetic, potential, and interaction energy. The 
interaction is split into two terms, one H^^^^^^ only involves particles with the same crystal 
momentum, while H^^^ involves particles with different crystal momentum. In (P) through 
(1) the argument k in ipj{k) is omitted. In @ the sum and integral are taken over all ki,ji 
such that momentum is conserved 

ki + k2-k3-k4 + 27r(ji + - js - J4) = 0, (10) 

and where not all of the ki are equal. In @ the shorthand notation ipi = ipj-lki) is used. 
The interaction is given hj g = g/2nL. 



The meaning of each of these terms is illustrated in FIG. 10. Solid lines show the kinetic 



energy of free particles as a function of the crystal momentum k. The periodic potential 
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conserves the crystal momentum and therefore only induces vertical transitions. The main 
effect of Hpot is therefore to split the degeneracies at the level crossings, giving rise to the 
band structure shown in gray. The two interaction terms scatter particles between these 
states. The 'vertical' interaction H^^^'^^\ only involves particles that all share the same 
crystal momentum, as illustrated on the right of FIG.pHjb. All other scattering processes, 
such as the one on the left of FIG.|ID|b are included in Hlj^. 

In our subsequent analysis we will ignore H[^. This is a quite drastic approximation 
which clearly restricts the phenomena which can be studied. For example, the superfluid- 
insulator transition seen in cannot be studied in this model. However, all studies of 
atoms in optical lattices which rely upon mean field theory (the Gross-Pitaevskii equation) 



implicitly make this truncation whenever they limit themselves to a Bloch ansatz This 
approximation therefore has a range of validity which is a superset of the mean field theory's. 
In particular this approximation is good when the interaction strength is the smallest energy 
in the problem. 

Once is eliminated, the sectors of different k are independent. If one identifies k with 
the $ in (^, then the two Hamiltonians are identical. For the remainder of the paper, we 
work with (^, while keeping in mind that all results can also be applied to a gas of particles 
in a periodic lattice. 



III. ENERGY LANDSCAPE IN ABSENCE OF IMPURITY 

Here, and in the next section, we calculate the properties of the microscopic model @, 
finding the general structures discussed in the introductory sections. We divide the discus- 
sion into several sections, based upon the limits of various parameters and the mathematical 
techniques used. 



A. two- mode approximation 

It is instructive to first analyze @ in the absence of an impurity, and in the limit where 
the interactions are sufficiently weak, ie. where A = and gN ^ 1, with N being the 
number of particles. The non-interacting single-particle energy states are shown as thick 
grey lines in FIG. ^ as a function of $. This spectrum and the physical properties which 
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we are interested in are periodic in $, and it suffices to consider — vr < $ < vr. The ground 
state, in the absence of interactions, consists of all particles condensed in the lowest energy 
state. Aside from providing a global shift in the chemical potential, weak interactions only 
introduce a significant perturbation when the energy difference between two levels is less 
than gN/2, so that the interactions mix the two states. Focusing on the level crossing at 
$ = vr, the system is reduced to two-levels with an effective Hamiltonian 



H^ = {<P + T^friQ + (0 - Tifm + + nl + 4noni), (11) 

where = $ — tt, and Uj = CjCj are constants of motion. The eigenstates \nQ,ni) = 
(co)"°(c|)"^ |0) /y/riolriil, have a fixed number of particles in each momentum state. Just 
like in the noninteracting system, for < (respectively > 0) the ground state is |A^, 0) 
(|0,A^)). Interactions play a role here only through the fact that when 7 = 2TT\(f)\/g < 
N/2 a barrier, as illustrated in FIG. [l^, separates these two states. This barrier exists 
because density modulations are required to transfer particles between angular momentum 
states. In the presence of interactions these modulations cost energy. The barrier state, 
\N/2 + 7, N/2 - 7), has energy Ei, = {(p + 7r)\N/2 + 7) + (0 - 7i)\N/2 - 7) + SgN^A - gj^, 
compared with E±. = (0 ± 7r)^A^ + gN'^/2 for the other extrema. The maximum barrier 
height (occurring at = 0) is 6E = gN^/A. 

This same scenario is repeated at all other level crossings in FIG. |TT|. Thus, in this weakly 
interacting limit, one can take the eigenstates of EQ. (^ to be the 'Fock' states 



|no,ni,n_i---) = n4=|0), (12) 

where the occupation numbers rij obey the constraints rij > and J2j % = For large 
numbers of particles there is no approximation involved in thinking of the rij as continuous 
variables. 

We have already detailed the energy landscape when we truncate this space to two n^'s, 
and the space of allowed states consists of a line (the x-axis of the plot in FIG. |T2|). When 
three n^'s are included, the space is a triangle, and with four n^-'s it is a tetrahedron. The 
d-dimensional generalization of a triangle/tetrahedron is often called a simplex or a hyper- 
triangle, and within our approximations, the eigenstates of the Hamiltonian form an infinite 
dimensional simplex. The corners of this simplex are cusps in the energy topography. Cusps 
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play a role similar to saddle points, as they are each classified by the number of independent 
directions in which the energy decreases. 



As explicitly shown in FIG. |T2|for the case of two levels, there exists a range of $ for which 
one can find extrema in addition to these cusps. In FIG. the extrema are marked by the 
number of 'downward' directions. Loops are clearly visible around points where the non- 
interacting states cross. These loops get smaller when the interaction strength decreases. 
Taking g and $ as the control parameters, one reproduces the structure in FIG. ^. 

As in our generic discussion, these loops lead to hysteresis. Suppose our model 1-D gas 
starts in the state |A^, 0) at $ = 0, and is then accelerated so that $ is slightly larger than vr. 
The barrier will then keep the system from jumping into the new ground state |0, A^). Thus 
the presence of this barrier implies the existence of persistent currents. Of course, since we 
have considered only very small interactions, these currents only exist near $ = ±7r. 

B. Mean field theory 

A useful tool to further illucidate our model (|^) is mean field theory. Here we use a mean- 
field theory to show that the general structure of metastability and superfiuidity found within 
the two-mode approximation continues to be valid for larger interaction strengths. 

This discuss reveals several important points. (1) Even though a one-dimensional Bose 
gas is a Luttinger liquid and is usually not studied using mean field theory, we show below 
that mean field theory correctly describes the behavior of a one- dimensional Bose gas for a 
significant parameter range. The exact details of this parameter range is discussed below. (2) 
The extrema of the mean-field Hamiltonian are in one-to-one correspondence with the energy 
extrema in the many-body Hilbert space discussed in section [III A| . (3) In the regime where 
both the two-mode approximation and mean-field theory are applicable, the Bogoliubov 
excitation spectrum of the mean-field coincides with the exact low energy excitations of the 
many-body problem. This correspondence is well-known from the Bethe ansatz analysis of 
the one-dimensional Bose gas |T^ . 
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1. The two-mode regime 



We begin by considering mean field theory in the regime where the two-mode approxi- 
mation is vahd. For weakly interacting Bosons, mean field theory can be formulated as a 
variational method in which one assumes all of the particles are in the same single particle 
state. In the effective two-level Hilbert space described by the Hamiltonian in EQ. (|ll|), we 
consider wavefunctions of the form \a,j3) = (aaj + /3a|)^|0), where the variational param- 
eters a and f3 satisfy |ap -|- = 1, and |0) is the vacuum state containing no particles. 
For < (0 > 0), the energy is minimized by |a;p = 1, = (jap = 0, = 1), which 
(within the two-mode approximation) is the exact ground state as shown in SEC. p.11 A . 
When 101 < g{N — l)/47r, one also finds a mean-field state which is a maximum of the en- 
ergy, corresponding to the barrier state found previously. The topography of the mean-field 
energy landscape mirrors that of the exact eigenstates, and the barriers found within the 
mean field theory do not differ significantly from the exact barriers. 

It is quite striking that the mean field theory compares so favorably with the exact two- 
mode calculation, considering that the exact barrier state contains two, rather than one, 



condensates and is therefore referred to as "fragmented" |Q . The connection between these 
states is understood by noting that that the mean-field barrier corresponds to a 'phase slip', 
where the fiuid density vanishes at some point. By its nature such an event breaks rotational 
symmetry. Averaging over the possible location of the slip restores the broken symmetry, 
and leads to the exact (fragmented) barrier state ^3 . 



In this weakly interacting limit the excitations of the mean field theory correspond to 
the exact low-lying excited states of the two-mode system. This result is trivially obtained 
by substituting our variational wavefunction into (0) and calculating the frequencies of 
small oscillations. For those familiar with dilute gases of bosons, this result is perhaps even 
simpler to derive by going beyond the two mode approximation, and directly writing down 
the excitations of a condensate moving at velocity Vc- As with a three-dimensional system, 
an excitation of wave-vector k is given by the Bogoliubov form. 



Ek = — VA;2 (fc2 + 2gN/L-') - v,k (13) 

Here the finite size of the ring restricts the wavevector to A; = 27Tn/L, with integer n. 
Similarly Vc = 27ihn'/ {mL) is quantized with integer n' . In the frame rotating with velocity 
VL the excitations have energies E}^ — hQL /2tx. Within the two- mode approximation we are 
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limited to n = 1, and it is straightforward to verify that excitation spectrum matches the low 
energy spectrum calculated directly from (pA]). Exploring the excitation spectrum around 
the saddle-point state, one finds a zero-mode corresponding to translations of the phase slip, 
and negative energy modes corresponding to falling towards one of the local minima. In the 
exact two-mode theory, the zero mode corresponds to changing the relative phase between 
the modes. 



2. Beyond two-modes 

In addition to the insights provided above, the mean-field approach also provides a sys- 
tematic way to explore EQ. for interaction strengths which are beyond the scope of the 
two mode approximation. The mean-field theory involves replacing the field operators 
with c-numbers. It is convenient to work in real space, defining a 'condensate wavefunction' 

ip{x) by 

Ck = f^dxe'^^tPix), (14) 
Jo 

where / dx|^/'(x)p = A^, is the number of particles. An energy landscape can be found in 
the space of all possible square integrable complex- valued functions ip^x). This landscape 
contains all of the structures seen in the two-mode version of mean field theory. In particular 
the k = n states given by ipn{x) = v/iVe*^*^^, are always stationary points. Their stability 
is given by the Landau criterion that if the excitation spectrum in EQ. (|T3|) is positive 
then they are local minima. Otherwise there exists a direction of negative curvature. The 
existence of multiple minima in the energy landscape leads to hysteresis and superfluidity. 

When both the k = n and k = n — 1 states are locally stable, there must exists a 
saddle point separating them. As in the two-mode case this saddle point involves a 'phase 
slip', where the density vanishes and the number of units of circulation can change. The 
real space wavefunction of the phase slip takes on the form of a hyperbolic trigonometric 
function, whose exact form was determined by Langer and Ambegaokar in the context of 
superconductors The barrier height is understood by recognizing that the length-scale 



for the phase slip is the healing length ^ l/^/gN in our dimensionless units). The 
presence of the slip increases the density from A^ to N/ (1 — ^), at an energy cost per particle 



of gN^ h J agL/ d\/mL'^ in physical units (assuming transverse harmonic confinement 



with lengthscale d±_ = ^Ti/muj). A more careful calculation, detailed in |]2S|, verifies this 
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result with a coefficient w32/9. 



3. limits of validity 

We have shown that within mean-field theory the one-dimensional Bose gas is superfiuid in 
that it exhibits hysteresis under changing the rotation speed. It is therefore very important to 
understand the limits of validity of mean-field theory. We estimate these limits by calculating 



the depletion within a Bogoliubov approach (for example see ||26|), where one finds that 



where = ^^/2m(fc^ + gN/L'^). Only when this ratio is small compared to 1 is the mean- 
field theory valid. For experimentally relevant parameters, the logarithm is of order ten, and 



the prefactor y g/N determines the size of the depletion. A physical interpretation of this 
factor is that the healing length ^ {1/ y/gN) which governs the scale of the phase slip must 
be larger than the interparticle spacing (1/A^) for fiuctuations to be small and for mean-field 
theory to be applicable (ie. phase slips cost very little energy if they fit between particles). 
To study stronger interactions one needs to include short-range fiuctuations - either through 
'bosonizing' the system , or by using the Bethe ansatz . 



Substituting plausible experimental values into EQ. |T3], show that it is much easier to be 
in the regime where mean field theory is applicable than it is to be in the strongly correlated 
regime. For example, with 10^ atoms of ®^Rb (with scattering length as ~ 5nm) in a ring 
of circumference L = lOO/xm, and transverse confinement u!± = 500s^^ (corresponding to 
d± ~ lyum) one finds 6N/N ~ 0.2%. Decreasing N or increasing L leads to proportionally 
more depletion and can bring one into the strongly correlated regime. 



IV. PERSISTENT CURRENTS IN PRESENCE OF IMPURITY 

Having established the superfiuid behavior of a weakly interacting gas through analysis of 
the energy landscape, we now analyze the behavior of such a system when a small impurity 
is added. We shall see that as long as the impurity strength is small compared to the 
interactions, such an impurity leads to extremely slow (typically exponentially slow) decay 
of persistent currents. In the opposite limit, where the impurity is strong, no persistent 
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currents nor hysteresis exists. 

Even when the impurity is weak, it does have a dramatic effect on the (mean-field) energy 
landscape in FIG. |Tl|, in that it opens up gaps as seen in FIG. |ll]b)- In the case of atoms in a 
periodic potential these gaps are the familiar band-gaps from solid-state physics discussed in 
the introduction. Note that the gaps do not change the fact that one has hysteretic behavior 
signaled by the swallow-tail loops. In this section we present the quantitative theory of these 
loops in the limit of weak interactions and calculate the quantum tunneling from one local 
minimum to another (ie. the decay of persistent currents induced by quantum fluctuations). 



A. two-mode model 



Interactions are most easily understood within a two-mode model. For its validity we will 
need both weak interactions gN <^ 1, and a weak impurity A <^ 1. The ratio X/gN which 
compares the impurity strength to the interaction strength is arbitrary. As before, it suffices 
to consider the system near = $ — vr = 0, where the Hamiltonian may be truncated to 
(|ll|), with an additional impurity term i^imp = A(a|ao + ajai). This system of bosons in 
two states can be mapped onto the precession of a spin N/2 object obeying a Hamiltonian 

iJspin = 47r0^, + 2A^, - 2gSl + Cn (16) 



where 5*2 = {aQao — a\ai)/2, = {S+ + S-) /2, Sy = {S+ — S-)/2i, S+ = alai, and 5*- = a\ao 
obey the standard spin algebra, and C^v is an uninteresting c-number. This mapping is the 
inverse of the method of "Schwinger bosons" . The quantum dynamics of such large spins 



are well understood |£8|, so we only briefly outline the analysis necessary to understand the 



mean-field structure, and the decay of persistent currents. 



The classical (mean field) energy landscape of the spin is shown in FIG. |T3|, where the 
angles 6 and describe the direction in which the 'spin' is pointing. In the absence of 
interactions, there are two stationary points, a maximum and a minimum. These represent 
the first and second band in the single-particle energy spectrum. For sufficiently strong 
interactions 

{gNf/' > (27r0)2/3 + A^/^, (17) 

a second minimum and a saddle spontaneously appear (a saddle-node bifurcation), giving 
rise to the swallow-tail energy spectrum shown in FIG. [Tl|. The new minimum is analogous 
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to the metastable states found in the absence of the impurity. The difference here is that this 
metastable state is distinct from the upper band, which is still represented by the maximum. 
In the symmetric case (0 = 0), the minima occur at = — (X/gNy, and the barrier 

has a height Eb = gN'^/2{l — X/gN)"^ as is clear from (p!^), the metastable state can only 
exist if the interaction strength gN exceeds the impurity strength A. 



B. Quantum tunneling in the two- mode limit 

Quantum mechanically, there are matrix elements for tunneling from the upper minima 
to the lower, and the upper state acquires a finite lifetime. For those more familiar with 
particle tunneling than with spin tunneling it may be helpful to instead map the problem onto 
the motion of a particle on a one-dimensional lattice. Introducing states \m) = (— l)™'|r;,o = 
N/2+m, rii = N/2—m) and operators c\^, which create these states, the many-body problem 
in the two-mode approximation is equivalent to a single particle with a Hamiltonian 



m - gm-) c^^c^ 



N 



AW ( Y + + l) - {cln+lCm + cl,Cm+l) 



(18) 

where, C'^ is an uninteresting c-number. This Hamiltonian represents a particle on a lattice 
with an inverted parabolic potential and an unusual spatially dependent hopping (which can 
be viewed as a spatially dependent mass). A particle surmounting a barrier by hopping on 
such a discrete lattice behaves somewhat differently than a similar particle with a continuous 
coordinate. In particular, for sufficiently large barriers, the tunneling rate is a power-law in 
A rather than the familiar exponential. Such power-law behavior was predicted by Kagan 
et al. in discussing the finite temperature lifetime of persistent currents . 

An important conceptual point to consider here is to what extent the scenario discussed 
so far can lead to the decay of a current. We have reduced the many-body problem to 
the one-dimensional quantum mechanical motion of a single particle. There is no source 
of dissipation within this model, and one would naively expect to see coherent oscillations 
between the two wells rather than a decay. This intuition is in fact correct when the two 
wells are degenerate (ie. = 0), or when the tunneling is extremely weak. This coherent 



tunneling limit is illustrated in FIG. |lja, in which only the lowest state in each well is 



19 



coupled. When there is a significant mismatch in the energies of the two wells, say the 
left hand well has more energy than the other, the situation is different. As pictured in 
FIG. the ground state on the left can be coupled to several excited states on the right, 
a situation analogous to an excited atom coupled to a large number of vacuum modes. In 
the limit of large A^, the spacing between the modes in the right hand well vanish, and the 
state on the left hand side of the barrier acquires a finite lifetime. A detailed study of this 



cross-over from coherent oscillations to decay is found in [^. One expects decay whenever 
the tunneling rate is significantly larger than A/^, where A is the characteristic energy 
spacing (here given by the inverse of the density of states of Bogoliubov phonons). 

A semiclassical analysis of (IS) is detailed in [^|. The eigenvalues of (0) are found by 
solving a difference equation. This difference equation can be approximated by a differential 
equation which is amenable to a WKB analysis. The resulting expression for the lifetime r 
of the persistent current is 

r = roexp ( / arccosh — ^ , ^ ds| (19) 

\Jbi 2\J{N/2f - s2 ; 



The attempt frequency t^^ is roughly the frequency of small oscillations about the local 
minimum of ([T6|). This frequency also coincides with the quantum mechanical energy of 
excitations in the metastable state (which as already pointed out, is the energy of Bogoliubov 



excitations). For the symmetric case, = 0, one finds Tq"^ ~ 2y {gNY — A^. The integration 
variable s corresponds to the projection of spin Sz-, the limits of integration hi and 62 are the 
classical turning points, and E is the energy, given by (p^ with the constant Cn removed. 
When s = 61^2 the argument of the arccosh is 1. 

For sufficiently small barriers, ^ AA^, the arccosh can be expanded as arccosh(l + 
x) = \/2x + 0{x^^^), yielding a decay rate which is exponentially small in the impurity 



strength, r ^ ~ e'xp{—ayNEh/X), where a is of order unity. For larger barriers, one 
can expand arccosh(z) = log(22;) + 0{z~'^), which leads to a power law behavior, 
{NX/ Eb)~^^''^^ , where /5 ~ 1. Unless the barrier is tuned extremely close to zero, both 
of these expressions yield astronomically large lifetimes whenever a barrier exists. Thus, 



the condition for superfluidity reduces to the condition for a barrier (|T^. In FIG. |T^, 
equation (0) is numerically integrated for some representative parameters, verifying these 
asymptotics. 
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FIG. 1: Typical energy landscapes: energy E vs order parameter M. In a) two minima (labeled 
1 and 2) are separated by a barrier (3). In b) one minimum and the barrier coalesce. In c) only 
one minima exists. A control parameter (H) tunes from one landscape to another. A plot of the 
energy of the extrema versus the control parameter is shown in FIG. 




FIG. 2: Energy extrema as a function of control parameter H. Solid lines denote minima, dotted 
denote maxima/saddles. The points labeled a, b, and c coincide to the energy landscapes in figure 
which respectively have 3, 2, and 1 extrema. The points labeled 1, 2, and 3, coincide with the same 
points in FIG. |l]a. The existence of multiple minima at the same value of the control parameter 
is a ubiquitous sign of hysteresis. The presence of two minima requires a maximum/saddle (see 
FIG. HI). 
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FIG. 3: Three dimensional depiction of a swallow-tail energy structure. The x and y axes represent 
control parameters, here referred to as Hi and H2, while the vertical axis is the energy E. The 
self-intersecting surface shows the stationary points of the energy. A two-dimensional slice is seen 
in FIG. where Hi can be identified with H. 
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FIG. 4: Catastrophe set: values of control parameters for which the number of extrema of the 
energy structure in FIG. ^ change. Inside the cusp there are three extrema (two minima and a 
maximum) while outside there is only one. 
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FIG. 5: Hysteresis loop in a superfluid. The rotation rate of the fluid is shown as a function of 
the rotation of the container. The arrows denote the direction of the hysteresis loop. As seen in 
experiments on helium the sloped lines are actually made up of many discrete jumps, which 
cannot be resolved on this scale. 
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FIG. 6: Band structure. The solid line shows the energy of Bloch waves of crystal momentum k 
for non-interacting particles in a periodic potential within an extended zone scheme. The dashed 
lines show the energy states in the absence of the periodic potential. 
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FIG. 7: Schematic of energy extrema for a condensate in an optical lattice. Solid and dotted 
lines denote minima and saddle points. One of the spinodal points, where the number of extrema 
change, is marked by a grey circle. 




FIG. 8: Geometry of a double-well trap. The potential energy V is shown as a function of a spatial 
coordinate r. The two wells have an energy difference 5. 



26 




FIG. 9: Illustration of hysteresis in a double well trap filled with attractive bosons. In each picture 
the trap from FIG. ^ is shown with a set value of 6, and the particles shown by small grey circles 
in one of the two wells. In a) the bias is positive, and all of the particles are in the left well. The 
bias is slowly switched to a small negative value in b). The particles remain in the left hand well, 
even though the ground state has all of them on the right. For large enough negative detuning c), 
the particles all jump to the right. 
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FIG. 10: Illustration of the terms in Hamiltonian (^. In a), the black line shows the kinetic 
energy i/kin of single particle states as a function of crystal momentum k. The periodic potential 
Hpot couples states with the same /c, splitting the degeneracies, giving rise to the band structure 
shown as a dotted grey line. In b) scattering processes are illustrated. On the left hand side 
a generic scattering event is shown where two particles with arbitrary momenta, scatter to two 
other states. On the right a 'vertical' scattering event is shown, where two particles with the same 
crystal momentum scatter to two other states, preserving k. These vertical scatterings are included 
in hI^^^'^ , while all others are in H[^^ 
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FIG. 11: Spectra: Thick grey lines lines - single particle energy levels of non-interacting particles 
in a 1-D ring of length L. Black lines - mean-field energy extrema of interacting system, with 
global mean-field shift removed. Solid are local minima, dotted (dashed) are saddle points/cusps 
with one (two) direction of negative curvature/slope. Figures a) and b) are respectively with and 
without an added impurity. In b) the area around <^ = vr is enlarged and displayed in an inset. 
Notice the similarities between the single particle states for particles in the ring shown here, and 
the band structure for particles in a periodic potential in FIG. |l^ 
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FIG. 12: Energy barrier separating persistent current carrying states within the two-mode ap- 
proximation. The ordinate shows no, the number of particles in the / = state, the remaining 
N — uq particles are in the / = 1 state. The maximum occurs at uq = uq = N/2 + 2Tr(j)/g, where 
is a measure of the rotation speed, and g the interaction strength. A barrier only exists when 
< no > A^. The labeled energies Ei, and E± are given in the text. 
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FIG. 13: Energy landscape of spin which represents the state of the one dimensional Bose gas 
within a two-mode approximation in the presence of an impurity. The number of particles with 
angular momentum / = (/ = ?i) is A^cos^(^) {Nsm'^{9)), while 4> is the phase angle describing 
the coherence between these states. For the contour plots, a) and b), darker colors represent lower 
energy; stars, triangles, diamonds, and squares represent the global minimum, global maximum, 
local minimum and saddles. In a) the inequality in EQ. ( [l7| ) is not satisfied, and only two extrema 
exist, while in b) the inequality is satisfied. A nonlinear scale is used for the contours in b) to 
emphasize the extrema. All of the extrema occur on a great circle parameterized by setting = 
and letting 9 run from to 27r. This corresponds to a path from the 'south pole' to the 'north 
pole' along the 'meridian' where = 0, and returning along the = vr. Figures c) and d) show the 
energies of a) and b) as a function of 9 along these paths. The structures from FIG. || are clearly 
seen. 
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FIG. 14: Schematic of tunneling within a double weh system: a) coherent tunneling between 
resonant states; b) incoherent decay from a single state on the left to a large number of states on 
the right. (Here this decay corresponds to a transition from a current carrying state to a stationary 
one with a large number of phonon excitations.) 
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FIG. 15: Lifetime of a metastable current carrying state (|T9|) in a toroidally trapped Bose gas 
in the presence of an impurity of strength A. For this figure, we consider the case where the trap 
is rotating at exactly one half quantum of circulation (p = 0. In the curves shown, the trap is 
described by a circumference L, a transverse confinement frequency uj, a number of particles A'^. 
These are respectively: solid line, L = lO^^m, N = 10^, to = 500s~^; dashed line, L = 10~^m, 
N = 10'^, UJ = 500s-^ dotted line, L = W'^m, N = W^, oj = SOOOOs^^ One expects that 
imperfections in the apparatus would lead to A ^ fiu). In all cases the particles are taken to have 
scattering length = 5nm and mass m =85 atomic units. The metastability turns on so quickly 
that the curves appear nearly vertical, even on this logarithmic scale. 
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